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Rapid Star Formation and Global Gravitational Collapse 
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1 INTRODUCTION 



Roughly a decad e ago, Ballesteros-Paredes et al.l |l999) and 
Hartmann et al. (2001) argued that star-forming molecu- 



lar clouds in the solar neighborhood evolve rapidly and 
produce stars on short - dynamical - timescales (see also 
lElmegreenl l2QO0h . The starting point for this picture was 
the observation that most nearby molecular clouds of sig- 
nificant mass are forming stars with typical ages of ~ 1 - 2 
Myr; only a small fraction of the stellar population exhibits 
ages > 5 — 10 Myr. A straightforward interpretation of the 
observations is that local star formation ensues quite quickly 
after molecular cloud formation, and that lifetimes of these 
nearby star-forming clouds are typically only a few Myr. 
Furthermore, in some cases the spread in ages of the bulk 
of the stellar population was considerably less than a lateral 
crossing time. To explain these observations, we proposed 
that molecular clouds in the solar neighborhood tend to be 
formed by "large-scale flows", accumulating material in a 
direction r oughly perpendicular to the lateral extension of 
the c l oud (iBallesteros-Paredes et al] 1999; H artmann et al.l 
l200ll : iHeitsch et alJl2008bh . Building clouds in this manner 
thus eliminates the need for communicating the "informa- 
tion" needed to trigger star formation roughly simultane- 



ABSTRACT 

Most young stars in nearby molecular clouds have estimated ages of 1-2 Myr, suggest- 
ing that star formation is rapid. However, small numbers of stars in these regions with 
inferred ages of > 5 — 10 Myr have been cited to argue that star formation is instead a 
slow, quasi-static process. When considering these alternative pictures it is important 
to recognize that the age spread in a given star-forming cloud is necessarily an upper 
limit to the timescales of local collapse, as not all spatially-distinct regions will start 
contracting at precisely the same instant. Moreover, star- forming clouds may dynami- 
cally evolve on timescales of a few Myr; in particular, global gravitational contraction 
will tend to yield increasing star formation rates with time due to generally increasing 
local gas densities. We show that two different numerical simulations of dynamic, flow- 
driven molecular cloud formation and evolution 1) predict age spreads for the main 
stellar population roughly consistent with observations, and 2) raise the possibility 
of forming small numbers of stars early in cloud evolution, before global contraction 
concentrates the gas and the bulk of the stellar population is produced. In general, the 
existence of a small number of older stars among a generally much- younger population 
is consistent with the picture of dynamic star formation, and may even provide clues 
to the time evolution of star-forming clouds. 
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ously along the length of the cloud. The swept-up material 
is initially atomic; only after substantial column densities 
develop, as a result of both accumulation of gas and lateral 
gravi tational contractio n, does the cloud become molecu- 
lar ([Bergin et al.l [2004^ . This evolution, driven largely by 
gravity at late stages, helps explain why star formation is 
initiated shortly after "molecular cloud formation" Q 

Since then there has been substantial dis cussion of ap- 
parent age spreads in star- forming regions fsee lJeffriesll201ll : 
Jeffri es et al.l l201lh . One of the common findings is that, 
even though the bulk of the stellar population is young, 
there exists a small number of stars with apparent ages 
^5 — 10 Myr or mor e which seem to be mem bers of the re- 
gion (see for example IPalla et aL 1 120051 . l2007h . The question 
then arises: does the presence of a few older members in star- 
forming clouds contradict the idea of dynamic cloud evolu- 
tion and star formation? Do these apparently older stars 



1 HBB01 specifically limited their discussion to the solar neigh- 
borhood, where most of the gas is atomic and therefore molec- 
ular clouds must be made from atomic gas. In other regions, 
where most of the gas is molecular, there should be more non- 
star-forming molecular clouds. 
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Figure 1. Histograms of the estimated ages of stars in the ONC, 
from Da Rio et al. (2010), using isochrones from Siess et al. (2000) 
(left), and Palla & Stahler (2000, right). The observations, binned 
in linear rather than the typical log age, show a strong skewed 
behavior with time (see text) 



instead indicate a long phase of quasistatic cloud evolution, 
possibly supported by turbulence and/or magnetic fields? 

In this paper we use numerical simulations to show 
that dynamic models of cloud formation can account for 
age spreads comparable to those observed, without intro- 
ducing turbulent or magnetic field support, simply be- 
cause fluctuations in initial conditions as a function of po- 
sition result in some regions collapsing faster than oth- 
ers. Moreover, the dynamic models exhibit global gravita- 
tional collapse, which produces an increasing rate of pro- 
tostellar core and star formation over time, in a man- 
ner qualitatively similar to the acce lerating star formation 
rates inferred from observations by Pall a &; Stahlerl {l999, 
2000 ). A similar argument has been made i ndependently by 
Zam ora-Aviles &; Vazquez- Semadenil ([201 lh , on a semiana- 
lytical basis. 



2 A CASE STUDY; THE ORION NEBULA 
CLUSTER 

To illustrate the issues typically presented by observations 
of young stellar populations, we use results from recent ob- 
servational analyses of the Orion Nebula Cluster (ONC). We 
focus on the ONC because it has been so well- studied, and 
because its density allows one to limit consideration to a 
restricted area of the sky, thus minimizing the problems of 
possible contamination (see §4). 

Figure [J show s the age distribution determined by 
|Pa Rio et al.l (|2010h using two differing se ts of evolution- 
ary tr acks, and binning linearly in age as in IPalla fe Stahlel 
(2000). (We use only the observed members of Da Rio et 
al. without correction for completeness, but the difference 
is negligible for our argument.) Based on their analysis, 
the majori t y of the stars have ages < 4 Myr using the 
Siess et all d2QQ0h evol utionary tracks, or < 2 Myr using the 



Pall a fe Stahlerl dl999h t r acks (Figure [p. In a later study of 
the data, Reg giani et al.l ([201 lh infer that the ONC stars are 
not coeval, with star formation activit y between ~ 1 .5 an d 
3.5 Myr. Similar results were found by 1 Jeffries et al ] (|201lh . 
who estimated that the apparent mean age of the ONC is 
about 2.5 Myr, with 95% of the low- mass stars formed be- 
tween 1.3 and 4.8 Myr. 



It is important to consider these estimates in the con- 
text of the crossing time of the region. The Da Rio et al. 
observations span a region of about 30 arcmin north-south, 
or about 3.6 pc at the adopted distance of 414 pc. The 
(one dimensional) veloc ity dispersion of stars in the region 
is ab out 2.5-3kms _1 (j Jones fe Walker! 1 1 988l : iFuresz etaD 
2008), implying a full crossing time of ~ 1.2 — 1.4 Myr. 
The age spreads of the main peaks in the stellar distribu- 
tion (~ 4 Myr using the Siess et al. tracks; ~ 2 Myr using 
the Palla & Stahler tracks) are thus ^2 — 3 crossing times. 
These values do not create a major difficulty for the picture 
of rapid star formation, as the onset of gravitational collapse 
of individual objects need not be coordinated be tter than a 
small number of dynamical timescales (see, e.g., Elme greenl 
2000). Moreover, it should be emphasized that these dynam- 
ical timescales refer to the current state of the region. If, as 
we suggest in §4, the ONC region has contracted significantly 
over the last few Myr, the relevant dynamical timescale is 
longer than the crossing times estimated above. 

Finally, it is worth noting that observational uncer- 
tainties and problems with theoretical isochrones can pro- 
duce spu rious age spreads of similar orde r to those discussed 
above (|Hartmannll200lL l2003h . Indeed. [Jeffries et al.l (|201lh 
suggest that this apparent age spread is dominated by a 
combination of observational uncertainties and differences 
in the formation processes of individual stars 

The biggest challenge to dynamic models of star and 
cloud formation is the "tail" of o lder stars with appar- 
ent ages between - 5 and 10 Myr. Palla fe Stahlerl (|20Q0h 
showed that skewed apparent age distributions such as in 
Figure [1] were typical of nearby star- forming regions, and ar- 
gued that this was evidence for accelerating star formation 
over periods of ~ 10 Myr in molecular clouds. There are 
observational problems which can lead to spurious large age 
spreads (§4), but these may not account for all of the appar- 
ently older stars, especially objects with infrared excesses 
implying the presence of circumstellar disks, or signatures 
of accretion (jPalla et al.ll2007l : I Jeffries et alJfeOllh . Can the 
dynamic picture of star formation be reconciled with the 
presence of such older stars? 



3 STAR FORMATION WITH GLOBAL 
GRAVITATIONAL COLLAPSE 

In the dynamic picture, star-forming clouds are not in a 
quasi-steady state, but instead are continually evolving. 
Initially, the cloud forms by sweeping up mass via large- 
scale flows driven by stellar energy input or perhaps spi- 
ral density waves; eventually, gravitatio nal collapse leads to 



runaway contraction in local regions jHeitsch etaD I2006L 
^' " 120071: iHennebelle et all 



2008b; Vazquez- Semadeni et al. _ _ _ . , 
2008: 1rleitsch fe Hartmann! 120081 : iHeitsch et alJl2009h . T~ 
bulence leads to density fluctuations which are the seeds of 
subsequent gravitational collapse, modified by global con- 
traction which increases densities everywhere prior to dis- 
ruption by outflows, stellar winds, supernovae, etc. Because 
some initial fluctuations will create denser structures than 
others, some regions will collapse before others, even without 
special turbulent or magnetic support; and because cloud 
densities increase with time, initially due to accumulation 
of material in the post-shock regions, and later due to grav- 
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Figure 2. History of massive core formation in the simulation 
Gs of Heitsch et al. (2008b), plotted as a function of time prior 
to the end of the simulation (see text) 



itational collapse, one would expect the star formation rate 
to increase with time. 

To illustrate this sequence, we examine the results of 
two numerical simulations of cloud formation from colliding 
large-scale atomic flows. (The mechanism is more general 
than the specific setup designed for computational conve- 
nience, as any swept-up flow can be turned into a colliding 
flow in a frame of reference moving with the swept-up gas.) 
It should be emphasized that these simulations minimize dif- 
ferences in the onset of collapse along the cloud because the 
flows collide all along the interface exactly at the same time. 
Moreover, no injected turbulent support or magnetic fields 
are included. Even in these idealized cases, star formation 
occurs over a finite time. 

The firs t simu lations we consider are those of 
Heitsch et al. | (|2008bh . At a resolution of ~ 0.08 pc, these 
calculations cannot follow fragmentation down to (low- 
mass) protostars; instead, as a proxy for star formation, 
we use the criterion of the formation of gravitationally- 
bound, dense cores. W e identify cores initially via clumpf ind 
([Williams et al.ll 19941 ). We consider only gas with T < 100 K 
to avoid having to analyze all the volume of the simulation, 
and demand that the local free-fall time in the clump is at 
least a factor of ten smaller than the global free-fall time of 
the cloud. The bulk of the core masses have temperatures 
near 20 K. 

Figure [2] shows the ti me sequence o f massiv e core for- 
mation in simulation Gs of IHeitsch et al.l (|2008bh . As these 
cores are quite massive, we alternatively plot both the num- 
ber of cores and the mass in cores as a function of "look back 
time" from the end of the simulation. Although the reso- 
lution is limited and the timesteps are crudely binned for 
reasons of statistics, it is clear that this model produces an 
accelerating rate of core formation (as well as mass growth). 

Figure [3] shows what is happening globally. After a suf- 




Figure 3. Surface density plots for model Gs (see Figure as a 
function of time, seen face-on (upper set of panels) and edge-on 
(lower set). The colors correspond to the column density in log 
cm 2 as given in the color bars. 



ficient amount of mass is accumulated in the post-shock re- 
gion by the colliding flows, gravity begins to take over, re- 
sulting in a more rapid increase in densities. The structures 
are m ore local than in sma ller cloud simulations Gfl and 
Gf2 of IHeitsch et al.l (|2008bl ). where the densest regions are 
the result of collapse into a filament. The latter simulations, 
though harder to use for tracking core formation because 
less mass was involved and therefore fewer cores were made, 
illustrate the type of evolution we envisage for the ONC, 
where gravitational contraction over a few Myr has lead to 
the dense gas residing in the narrow "integral-shaped fila- 
ment" (§4). 

The cores in this simulation are quite massive, rang- 
ing from about 150 Mq to about 8OOM0. To relate this to 
star formation requires an assumption that these cores do 
not remain inert, but instead continue to collapse beyond 
what we can determine given our resolution, fragmenting 
into more typical prot ostellar cores of a few solar masses 
(e.g. Andre et alJ feoiO) and then into stars. If the efficiency 
of star formation is directly related to th e amount of dense 
core mass, as seems reasonable (see, e.g. lLada et al1l2010l . 
and references therein) , then the core mass and number evo- 
lution of model Gs can serve as a proxy for star formation. 

The other simulat i ons w e consider are those by 
IVazquez-Semadeni et al.l ((2007), which are comparable to 
the previous ones, in that two warm, thermally bistable 
streams collide to form a dense, collapsing, cold cloud. These 
simulations used gadget, the s mooth Particle Hydrodynam- 
ics (SPH) code developed by Spring el et al . (2001), which 
includes the possibility of sink particle formation. 

In Figure |4] we show the back-in-time histogram 
of the newborn sink particle s for the fiducial run 20 
([Vazquez-Semadeni et al J l2QQ7h . In this figure, t = Myr 
corresponds to t ~ 23 Myr since the begining of the colli- 
sion and thus t — 6 Myr corresponds to the time in which 
the first sink particles are formed, i.e., to ~18 Myr from the 
start of the simulation(3 From Figure [4] it is clear that, as 
well as in the case of the dense cores in the previous figure, 



2 It is worth emphasizing that during the initial evolutionary 
stages the flow-formed cloud is atomic, and will only become a 
molecular (CO) cloud when column densities become sufficiently 
large to shield the molecules from the diss ociating interstellar 
radiation field; see (see Heitsc h fc Hartmannl uOOSV 
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Figure 4. Histograms of the newborn sink particle s as a func- 
tion of time for run 20 in Vazquez-Semadeni et al . (20071). The 
number of newborn sink particles increases rapidly with time 
during the first 4-6 Myr. After this time, the mass involved in 
sink particles is such that, assuming a typical IMF, there would 
be enough massive stars to d isrupt the cloud ((see for details 
IVazqu ez-Semad eni et al.1l2007l) 



the sink particles are being formed at an accelerated rate. 
We must note that at the final epoch plotted, the mass in- 
volved in sink particles is such th at there would be enough 
OB stars to disrupt the cloud (see IVazquez-Semadeni et all 
120071 ). assuming a typical IMF (e.g~ lKroupall200lh . Thus 
we omit the subsequent evolution of the simulations as likely 
being unrealistic. 

In both cases, the formation of dense structures (cores, 
sinks) occurs at an accelerated rate due to gravitational col- 
lapse, becoming denser and t hus evolving loca lly faster, es- 
pecially in filaments (see, e.g. JPon et al.|[201lh . While there 
are limitations to these simulations - for example, there 
probably would be continuing fragmentation beyond the res- 
olution limits, and the initial conditions are idealized - the 
qualitative resemblance of the simulations to the observa- 
tional results is suggestive. 



4 DISCUSSION 

Given the current state of observational constraints, pos- 
sible errors, and uncertain physics of star formation, one 
must exercise caution i n interpreting the age distributions 
of star-forming regi ons (|Hartmannll200lL 120031: 1 JeffriedlioTll ; 
Ijeffries et aP boilh . We agree with Ijeffries et aP (|201lh 
that it is especially difficult to be certain of age spreads 
in the range of ~ 1 — 3 Myr, particularly as differing initial 
conditions for protostellar formation may result in signif- 
icant initial dispersio n s in luminosity fearaffe et al. 2009; 
iHosokawa et alj 1201 ll ; lHartmann et al] 201lh . There may 
even be difficulties in assigning the ages of ~ 10 Myr-old 
stars; contamination by foreground stars is an issue, espe- 
cially for stars without disks or accretion, as star-forming 
regions are spatially-correlated. Even stars with disks can 
appear anomalously faint for their colors if observed edge- 
on, and thus mostly detected in scattered light. 

Nevertheless, our simulations show that the presence 
of small numbers of older stellar members in molecular 
clouds does not pose a particular problem for the idea of 
rapid or dynamic star formation. All such models begin 



with turbulent fluctuations, and it is plausi ble that a few 
espec ially dense perturbations collapse first (Hei tsch et all 
2008b). The effects of global gravity then generally results 
in ever- increasing densitie s, with runaway co n tracti on in 
su bregions, as argued by iBurkert &; Hartmann! (|2004T ) , and 
bv lHartmann fe Burkertl <|2Q07h specifically for the Orion A 
complex. In this scenario, a small number of stars are formed 
by a few especially dense initial turbulent fluctuations before 
the overall collapse leads to the main phase of star forma- 
tion. The two simulations presented here suggest that time 
spans of 5 to 10 Myr can be accomodated by purely dynamic 
models, as long as the initial star formation rate is quite low 
(e.g., compare Figures 2 and 4 with Figure 1). 

Absent observational problems, it is difficult to un- 
derstand the age distribution shown in Figure [1] with- 
out invoking substantial evolution of the ONC region over 
the last several My r . Fo r instance, if the suggestion of 
iKrumholz fe McKed ([2005) that the star formation rate per 
free-fall time is roughly a constant is correct, the age distri- 
butions in Figure [1] imply that the free-fall time has varied 
by an order of magnitude, and thus the average density by 
two orders of magnitu de, over the last 10 My r. Indeed, the 
recent simulations bv IKrumholz et aL | (|201lh of an "ONC- 
like cluster" show strong evolution over < 10 5 yr, while a 
somewhat longer cont raction tim e scale is ex hibited by the 
cluste r simulations of iBate et al. I (l2003h and iBonnell et al] 

(HoTil). 

The suggestion of significant cloud evolutio n is also con- 
sistent with kinematic stud ies of the ONC stars ([Tobin et al] 
2009; Proszkow et al ] l2009h . which suggest that both the gas 
and stars in the ONC are collapsing toward the central re- 
gions. More over, the spatial distrib ution of the stars (see, 
for instance, iDa Rio et aPl200l l2010h is wider in right as- 
cension than the narrow de nse "integral-shap ed filament" 
of molecular gas and dust (|Ballv et al] I"l987). This differ- 
ence is qualitatively consistent with global gravitational col- 
lapse; many stars could have been formed from the gas 
in the region in a more distended state, which has now 
collapsed to form a f ilame nt, as in the simple model of 
lHartmann fe Burkertl (|200/t ) for Orion A. 

The idea of large-scale gravitational collapse is also con- 
sistent with bo th observed column density pr obability den- 
sity functions (Ballest eros-Paredes et al] 201 lbh as well as 
with recent discussions showing that the "Larson laws" re- 
lating velocity dispersions ("turbulence ") with size scale s 
are not independent of surface density (|Hever et alj|2009h . 
These results can be interpreted as the natural outcome of 
star-forming molecular clouds bei ng in a state of hierarchical 
and ch aotic gravitational collapse (Balles teros-Paredes et all 
l2011ah . 

On the other hand, the concept of global collapse has 
been challenged by Dobbs et al. (2010), who argue on the 
basis of galactic-scale simulations that cloud-cloud collisions 
and stellar feedback prevent global gravitational forces from 
becoming dominant. In a narrow sense, this is not a problem 
for our picture, as we are focused on the dense star- forming 
regions of clouds, which Dobbs et al. agree do become bound 
(and form stars). Future observations of stellar proper mo- 
tions from the Gaia spacecraft might be able test whether 
or not star forming regions are globally collapsing. 
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5 CONCLUSIONS 

We have shown that models of flow-driven, dynamic dense 
cloud formation and evolution predict that star formation 
occurs over a finite interval of time necessarily greater 
than the timescales of local collapse for individual stars. 
Moreover, there is generally a strong increase in the dense 
core/star formation rate over time due to the increase in 
overall density and filament formation during global gravi- 
tational collapse. The small number of stars apparently older 
than a few Myr found in or projected upon star-forming re- 
gions may be a signature of this cloud evolution, though care 
must be taken to avoid observational problems. It may even- 
tually be possible to use carefully- vetted age distributions of 
pre-main sequence stars to infer the global evolution of the 
clouds from which they formed. 
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